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Introduction 


Advances in micro-fabrication processes have generated tremendous interests in miniaturizing 
chemical and biomedical analyses into integrated microsystems (Lab-on-Chip devices). To 
successfully design and operate the micro fluidics system, it is essential to understand the 
fundamental fluid flow phenomena when channel sizes are shrink to micron or even nano 
dimensions. One important phenomenon is the electro kinetic effect in micro/nano channels due 
to the existence of the electrical double layer (EDL) near a solid-liquid interface. Not only EDL 
is responsible for electro-osmosis pumping when an electric field parallel to the surface is 
imposed, EDL also causes extra flow resistance (the electro-viscous effect) and flow anomaly 
(such as early transition from laminar to turbulent flow) observed in pressure-driven 
microchannel flows. 

Modeling and simulation of electro-kinetic effects on micro flows poses significant numerical 
challenge due to the fact that the sizes of the double layer (10 nm up to microns) are very thin 
compared to channel width (can be up to 100’s of pm). Since the typical thickness of the double 
layer is extremely small compared to the channel width, it would be computationally very costly 
to capture the velocity profile inside the double layer by placing sufficient number of grid cells in 
the layer to resolve the velocity changes, especially in complex, 3-d geometries. Existing 
approaches using “slip” wall velocity and augmented double layer are difficult to use when the 
flow geometry is complicated, e.g. flow in a T-junction, X-junction, etc. In order to overcome 
the difficulties arising from those two approaches, we have developed a sub-grid integration 
method to properly account for the physics of the double layer. The integration approach can be 
used on simple or complicated flow geometries. Resolution of the double layer is not needed in 
this approach, and the effects of the double layer can be accounted for at the same time. With this 
approach, the numeric grid size can be much larger than the thickness of double layer. Presented 
in this report are a description of the approach, methodology for implementation and several 
validation simulations for micro flows. 


Description of the Sub-Grid Modeling 


To enable simulations without resolving the double layer, its effects on the flow can be modeled 
using the concept of a £ potential near the wall. A body force term, resulting from the interaction 
of the C, potential and the applied pressure drop, is added to the fluid momentum equations, 
written as: 
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The last term of equation (1) represents the electrokinetic effect. For streaming potential 
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( 2 ) 


induced by the applied pressure, 

K E i = ~ U iPe » 

and the local electric charge density, p e , is governed by the Poisson equation 

V 2 ^ = -^. (3) 

££ 0 

In the above equations, p is the fluid density, u; the Cartesian velocity components, p the 
pressure, Ty the stress tensor, s the fluid permittivity, 8=1 /k the Debye thickness of the double 
layer, and A 0 the bulk fluid conductivity. Under equilibrium situation, the ion distribution can be 

described by the exponential function, then equation (3) becomes the well-known Poisson- 
Boltzmann equation. For small surface charge situation investigated in this study, the linearized 
form of the Poisson-Boltzmann equation can be used. Equation (3) thus becomes: 

V 2 C = k 2 C (4) 

with BCs’: j^ 0 ’ y = 0 

K = °, y = °° 


C, is the potential field inside the electrical double layer and is the zeta potential. 

Equations (3) and (4) will be used to substituted into the last term of equation (1) for the 
electrokinetic effect. (We would combine £ and£’ 0 into one single symbol £ in the following 
formulation.) Considering the nature of C, potential, we analytically integrate the equation (4) 
combining with the boundary conditions, along a local wall normal direction to obtain the 
distribution: 

S = S 0 e“ Ky (5) 

Where, y is nonnal distance to the wall, and is the value of Q potential at the wall. The 
electrokinetic source term defined in equation (1) is only applied in a very narrow region near the 
wall. We analytically integrate it while neglecting the impact of C, potential outside the double 
layer, with assumption that the velocity profile within the double layer to be linear, and the 
application of Grahame equation [1] relating surface potential to surface charge density, <7, the 
source term becomes 


S = — a 2 fc-(l-3e~ 2 )Au (6) 

A 0 4 c 

Where, V is grid volume near the wall, A is the area of the wall face (see figure 1). The vector S 
denotes the momentum source contributions due to the streaming potential effect in different 
directions. The shear stress on the wall is calculated as: 
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t A Channel Wall 
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Figure 1. Integration Grid Volume. 

In deriving equation (7) we used the relation 5 = V . u c is cell center velocity of grid near the 

wall and u w is the wall velocity. The wall velocity is zero for a stationary wall. This definition of 
the effective distance essentially reduces the apparent wall shear stress, calculated based on the 
available information during calculations (i.e. u c , u w and 8). 

The above equations were implemented into ESI-Group’s CFD-ACE+ v2004 “Electro” module. 
Due to the sub-grid model nature, the grid system near the wall should be carefully constructed 
such that the first cell away from wall should cover the EDL thic kn ess. 

Results and Discussion 


The model was first compared to the benchmark test case of rectangular microchannel flow of 
Ren at al. [2], Detailed channel dimensions and testing conditions were described in [2], The 
model prediction and the experimental data are shown in Table 1. As discussed in [1], the 
electro-viscous effect would be significant when the ionic concentration of the solution is smaller 
the 5x10' M. As stated in [2], the data for KC1 aqueous solution of 10'“ M was taken to be the 
case without EDL effect. As can be seen from Table 1. The model predictions compared well 
(within 5%) with the experimental data. 



No EDL 

No EDL 

With EDL 
DIUF H20 

With EDL 
DIUF H20 

dP/dx (10 6 Pa/m) 

Re, pred. 

Re, exp. 

Re, pred. 

Re, exp. 

2 

1.01 

-1.05 

0.79 

-0.82 

3 

1.50 

-1.45 

1.20 

-1.25 

4 

1.98 

-1.90 

1.58 

-1.62 


Table 1. Prediction and experimental data [2] comparison. 


The second test case involved a T-junction micro flow on the Caliper N145 chip. The present 
model predicts up to about 40% deduction in flow rate when EDL effect is accounted for, as seen 
in Figure 2. Due to the uncertainties introduced by the dye in the experimental procedure for 
obtaining velocity data [3], the comparison is qualitative and the model shows correct trend for 
this complex micro fluidics system. 
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♦ KCI-2 Exp 
■ KCI-4 Exp 
KCI-2 pred 
KCI-4 pred 


Pressure Drop (psi) 


Figure 2. Prediction and data comparison for Caliper N 145 chip. 


Conclusion 


A sub-grid model was developed to account for electro-viscous effect in lab-on-chip 
microfluidics system. The model was successfully implemented into CFR-ACE+, and validation 
study was performed. The present model predicts straight micro channel flow well and 
qualitatively predicts T-junction flow correctly. The electro-viscous effects were found to be 
quite significant and reduce flow rate (up to 40% in T-junction flows), when compared to the 
classical laminar theory. Further refinement of the model to account for realistic velocity profile 
and non-linear Boltzmann distributions within EDL are desirable. Further validations with 
carefully designed complex micro flow experiments are also recommended. 
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